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Abstract. Construction of skeletonized path integrals for a particle moving 
on a curved spatial manifold is considered. As shown by DeWitt, Kuchaf 
and others, while the skeletonized configuration space action can be written 
unambiguously as a sum of Hamilton principal functions, different choices 
of the measure will lead to different Schrodinger equations. On the other 
hand, the Liouville measure provides a unique measure for a skeletonized 
phase space path integral, but there is a corresponding ambiguity in the 
skeletonization of a path through phase space. A family of skeletonization 
rules described by Kuchaf and referred to here as geodesic interpolation is 
discussed, and shown to behave poorly under the involution process, wherein 
intermediate points are removed by extremization of the skeletonized action. 
A new skeletonization rule, tangent interpolation, is defined and shown to 
possess the desired involution properties. 



1. SKELETONIZED PATH INTEGRALS 

The path integral method seems to provide a generic recipe for producing a quantum 
theory given a classical action. This recipe is based on the principle that amplitudes 
are given by sums, over the relevant histories, of the complex exponential of the action. 
For example, a wavefunction ilj{xi,ti) is propagated to a later time tf > ti by 

ip{xf,tf)= J V^e^'^^^/'^tPix^^U), (1) 
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where the integral is over all paths ending at the argument Xf. However, in order to 
provide a constructive definition for an expression like (|l|) one needs to give a meaning 
to formal concepts like an integral over paths, and it is in this step that decisions need 
to be made in the implementation of the quantum theory. 

The system considered in this work is a free non-relativistic particle of unit mass 
moving on a curved n-dimensional spatial manifold, which is described by the action 

5W = ^/rft,.,(x)^^. (2) 

This is the simplest action which exhibits the ambiguities considered here, which cor- 
respond in an operator theory to questions of operator ordering. One can also modify 
the action to include potential terms, but previous work [|l| suggests that the results 
will remain essentially the same. 

1.1. Configuration Space 

A formal path integral such as (|l|) can be given a concrete definition via the skeletoniza- 
tion process, in which the interval is broken up by a series of N time instants 
= 1, . . .iV}. The integral over all paths x{t) is realized, in the limit of an in- 
finitely fine time slicing, as a product of integrals over the positions {x(/)|J = 0, . . . N} 
of the particle at the discrete times: 




(3) 

where t(o) = U and t{Ar+i) = t/. This process replaces the action functional S[x\ with 
a function S/Lo 'S'(a;(/+i))f:(/+i)|x(/)t(7)) which is constructed by summing contributions 
from the intervals between time slices. Since the equations of motion are second order, 
there is a preferred path between a pair of endpoints (x, t) and (x',t'), namely, the 
classical path between them, which for the theory described by the action (|^) is an 
affinely parametrized geodesic.0 The action functional evaluated on that piece of the 
path is the Hamilton principal function S{x't'\xt). On the other hand, the measure 
ri/Lo/ ^{'^{i+i)t(i+i)\x(i)t(i))dPx(^i) for the path integral has no such natural definition 
in curved space. Feynman showed that A{x't'\xt) = [27iih{t' — t)]^"/^ was the 
appropriate measure for fiat space in Cartesian coordinates. The simplest generalization 
of the Feynman measure to curved space is 

Aoix't'lxt) = (4) 

where g is the determinant of the metric Qab at the point x, which ensures that the 
path integral is invariant under coordinate changes. However, there are other measures 
which are similarly invariant and also reduce to the Feynman measure when the spatial 
manifold is fiat. Expanding to second order in powers of the geodesic distance a' 
between x and x' (which turns out to be sufficient to evaluate the path integral), one 
family of such measures i^ 



^Throughout this paper, the points in question, even when separated by a finite distance, are taken 
to be close enough that there is no crossing of geodesies. 

^Specific measures of this sort were discussed in ||, ||], and placed into the family (||) in 



where y'^ are the Riemann normal coordinates of the point x' with respect to the origin 
X and i?"* is the Ricci tensor at x. 

This ambiguity in defining the path integral measure has tangible consequences. 
Namely, if we use a member of the family (|^) to define the measure in the skeletonized 
path integral (0), the wavefunction thus propagated satisfies 0, IT]] a Schrodinger equa- 
tion 

«^ = (-1 S"'V,V. + i^fi^fl) t). (6) 

The classical Hamiltonian g°'^PaPb is not quantized in the simplest way, by replacing it 
with an operator proportional to the covariant Laplacian, but also contains a measure- 
dependent term proportional to the scalar curvature. In particular, the straightforward 
generalization (Q) of the Feynman measure (in which A = 0) leads to a Schrodinger 
equation with a curvature correction term. 

1.2. Phase Space 

One way to resolve the measure ambiguity is to start instead with a phase space path 
integral 

ilj{xf,ts)= I D>P"a;e^^[P'^]/V(a;.,t.), (7) 
^/ 

in which the path integral measure is expressed as a product of Liouville measures 
dPp dT'x / {2Tih)'^ in the skeletonization 

= Jhn {\{J ^^^l^e^^(^(-^)*(-^)l^'^^^(^)*a.)/^) i^ix^^U). (8) 

Unfortunately, the ambiguity has now shifted into the definition of the contributions 
{S'(x(/+i)t(/+i) to the skeletonized action. This is because there is in general 

no classical trajectory through phase space satisfying the conditions 

x{t) = X, p{t) = p, x{t') = x', (9) 

since any two of these three conditions are sufficient to specify a classical phase space 
path. 

Since there is no classical phase space path, we need some interpolation rule which 
associates a virtual path through phase space with the set of boundary conditions (|]). 
(It should of course reproduce the classical path when the three boundary conditions 
are classically consistent.) A given interpolation rule defines a phase space principal 
function (PSPF) S{x't'\pxt) whose value for a particular set of arguments will be 
the action functional evaluated along the virtual path which the rule specifies for those 
boundary conditions. The interpolation rule then defines a skeletonization, under which 
and the skeletonized canonical action appearing in (P) is a sum of PSPFs. 

1.2.1. Geodesic Interpolation. One type of rule, introduced by Kuchaf pre- 
scribes the configuration space path to be the classical path from x to x', regardless 
of the value of p. The momentum then "goes along for the ride," being propagated 
in some linear way along the geodesic (Fig. |l]a). The resulting phase space principal 
function is ^ 

S{x't'\pxt) = pSa - ^-^G^\x'\x)paPb , (10) 




Figure 1. Geodesic interpolation (GI), illustrated in 1 + 1 dimensions, with the time direction 
running vertically and the space direction running horizontally, (a) Constructing a phase 
space path from the boundary conditions the configuration space trajectory is the classical 
path, a straight line from x to x', while the momentum p is held constant along the trajectory. 
In higher dimensions, if the spatial manifold is curved, the GI scheme also requires a rule for 
propagation of the momentum, (b) Failure of involution of paths for GI; using the values (|T9| ) 
of x' and p' which extremize the sum (18) of phase space principal functions, one constructs 



the two-step phase space trajectory, with the first step defined by x, p and x' and the second 
by x', p' and x" . Comparison to Fig. ||(a) shows that this is not the same path as would be 
constructed from x, p and x" by the GI prescription (although the two paths do have the 
same action when the space is flat). 

with the functional form of G°*(x'|x) depending on the details of the propagation rule. 
We call this "geodesic interpolation" (GI) because the property of the classical trajecto- 
ries which is retained by all of the virtual ones is that the configuration space projection 
of the path is a geodesic. This motivation for selecting rules in which the second order 
equation of motion, involving only the path x{t), is obeyed independent of the path 
p{t), comes from the invariance of the tensor formulation under point canonical trans- 
formations [x X{x),p P{p,x)], in which the transformation of the x components 
of a point in phase space is independent of its p components, but not under general 
canonical transformations [x — >■ X{p,x),p P{p,x)]. 

If the momentum p is taken to be parallel-propagated along the geodesic, the 
tensor (^"^(x'lx) in (|1^) becomes the inverse metric g""^ at x. This is one of a family of 
propagation rules discussed by Kuchaf which, for small intervals, have the expansion 

Gfix'lx) = g^' - ^R''^"'yWd + 0(a'=^), (11) 

where is the Riemann curvature tensor at x. 

Performing the quadratic integrals over the momenta in the skeletonized path 
integral (^) induces a measure 

A{x't'\xt)= , ^ (12) 

^J [27Tih{t' -t)Y' det G^^ 

on the remaining configuration space path integrals. The family of propagation rules 
given by ([TTD lead in this way to exactly the family of configuration space measures (|]). 
This means that they produce quantum theories whose Hamiltonians contain curvature 
correction terms ^^h^R; in particular, the parallel propagation rule (A = 0) leads to 



the Feynman measure (^) and hence to a curvature-corrected quantum Hamiltonian. 
The case A = 1, which Kuchaf showed arose from a propagation rule based on the 
equation of geodesic deviation, leads to a Hamiltonian with no curvature correction 
term. 

This work considers a criterion for selecting a phase space skeletonization rule, 
and hence a principal function. This criterion, as described in the following section, is 
based on a property of the PSPFs and paths themselves. 

2. THE INVOLUTION PROPERTY 

2.1. Configuration space 

We will judge phase space skeletonizations by whether they obey a property analo- 
gous to the involution property satisfied by configuration space skeletonizations. The 
configuration space property is defined as follows: Consider the contribution to the 
skeletonized action from the intervals between three consecutive time instants, called 
t, t' and t": 

S{x"t"\x't')+S{x't'\xt). (13) 

It is a function of three positions, x, x' and x", and has associated with it a path 
made up of two geodesies connected at the point x'. If we vary the location of that 
point x' while holding the endpoints x and x" fixed, the expression (|T3p is minimized 
by the x' which lies at the appropriate point on the geodesic from x to x". The path 
corresponding to that choice of x' is just the geodesic connecting the two endpoints x 
and x", and its action is the Hamilton principal function S{x"t"\xt): 

S{x"t"\x't') + S{x't'\xt) ^ S{x"t"\xt). (14) 

This property is known as involution. 

We say that the configuration space skeletonization procedure obeys involution of 
paths (lOP) because extremization of (|T^) with respect to the intermediate position x' 
leads to the classical path from x and x", and involution of functions (lOF) because 
([I^ ) holds. Clearly, the latter property is implied by the former. 

2.2. Phase space 

Given an interpolation rule, which defines a phase space principal function S{x't'\pxt), 
we can define the analogous property for a phase space skeletonization by considering 
two consecutive intervals in the skeletonization. The five quantities x, p, x', p' and x" 
determine a phase space path between times t and t" according to the interpolation 
rule. The action for this path is 

S{x"t"\p'x't') + S{x't'\pxt). (15) 

Extremization with respect to the intermediate variables x' and p' leaves a function of 
X, p and x"; we say that involution of functions holds if that is the same PSPF with 
which we started: 

S{x"t"\p'x't') + S{x't'\pxt) ^ S{x"t"\pxt). (16) 



The stricter condition of involution of paths is satisfied if the phase space path defined 
by the extremizing values of x' and p' is the same one which the interpolation rule 
would have constructed a priori from x, p and x" without making reference to the 
intermediate time slice t'. 

These concepts are easily illustrated in the case where the spatial manifold is 
flat; in that case the preferred propagation rule associated with geodesic interpolation 
is simply to keep the momentum constant (as measured in a Cartesian coordinate 
system). This leads to a PSPF 

t' -t 

S{x't'\pxt) = p ■ {x' — x) ^ — p^ . (17) 

The two-step action 

S{x"t"\p'x't') + S{x't'\pxt) =p' -x" + {p-p')-x' -p-x- ^—^p'^ - ^-^p" (18) 
is minimized by 

p = p and x' = x" — (t" — t')p. (19) 

Substituting these into (p!8|) gives p- (x" — a;) — ^^p^ = S{x"t"\pxt), so the G I scheme in 
flat-space does satisfy lOF. On the other hand, it is easy to see that it does not satisfy 
lOP, since the path obtained by extremization has a configuration space projection 
which starts at x at time t, follows a straight path to x" — (t" — t')p at t', and then 
makes a sharp turn before heading on to x" at t" (Fig. Qd). This is not the same as the 
path from x to x" which geodesic interpolation would dictate without the presence of 
the intermediate point, since that path is just a straight line from x to x" . 

If we look at the curved-space case in the limit that all of the distances are small, 
the zeroth-order results of course replicate the fiat-space ones (including the failure of 
lOP). If we consider the family of PSPFs defined by (pJ]), the behavior of the first 
correction terms will depend on the value of A. It turns out that for A = —1 (and for 
no other value) lOF continues to hold to the lowest non-trivial order. 

The reason why lOF can hold even when lOP fails is this: Since a given PSPF 
S{x't'\pxt) corresponds to the action functional evaluated along a virtual, rather than 
classical, phase space path satisfying the boundary conditions (||), there are many 
different paths which yield the same function. We can use this fact to our advantage 
by defining a different skeletonization procedure which produces the same PSPF over 
small intervals as A = — 1 geodesic interpolation, and thus also obeys the lOF property, 
but arises from a different family of paths, which satisfies lOP. 

3. TANGENT INTERPOLATION 

We will call the skeletonization procedure which respects lOP tangent interpolation 
(TI). As opposed to the geodesic interpolation schemes, which require the configuration 
space part of the path to be the classical geodesic, the configuration space path in the 
TI scheme is not even continuous. However, the classical equation relating momentum 
to velocity is respected by the prescribed path, so that the momentum vector remains 
tangent to the configuration space trajectory. Hence the name "tangent interpolation." 

The phase space path (Fig. ||a) determined by a given set of boundary conditions 
(P) begins at x with the momentum p at time t, jumps instantaneously to a point 




Figure 2. Tangent interpolation, illustrated in 1 + 1 dimensions, (a) The phase space path 
associated with the boundary conditions @ has an initial configuration space jump from x 
to X and then follows the classical phase space trajectory from x to x'; the point x is chosen 
so that the prescribed momentum p agrees with the initial momentum on the geodesic xx' . 
(In higher-dimensional curved space, this comparison is made in Riemann normal coordinates 
based at x.) (b) Involution of paths for TI. Using the x' and p' (|l9| ) which extremize the 
sum of phase space principal functions (p^), one constructs a two-step phase space trajectory, 
first from t to t' and then from t' to t" . There is no discontinuity at t' for this choice of x' 
and p' (i.e., x' = x'), and the initial jump from x to x implied by the chosen x' at time t' 
is the same as the one implied by x" at t" (along with the initial p in either case). Thus 
the two-step phase space path is the same as we would get by TI in one step from t to t". 
(Compare Fig. |2^.) 



X which is prescribed (in a manner detailed below) by the boundary conditions, and 
then follows the classical phase space trajectory from x at t to x' at t'. The point x is 
chosen so that in a Riemann normal coordinate (RNC) system based at x, there is no 
momentum discontinuity at the initial time t, i.e., so that the classical momentum at 
the beginning of the geodesic from x to x' has the same covariant RNC components as 
p. In flat space, this means that p = ^T^f , or 

X = x' -p{t' -t), (20) 

although in curved space the relation x{x't'\p x t) is only implicitly defined by the match- 
ing of momenta. However, the x which achieves the momentum matching condition can 
alternatively be found as the x which minimizes S{x't'\xpxt), the action along such a 
path, for a given x, p and x'. This action has two contributions, one from the initial 
discontinuity and one from geodesic from x to x', giving 

S{x't'\xpxt) = p^ya + S{x't'\xt). (21) 

Extremization with respect to x gives the curved-space counterpart of (|20|), and defines 
the PSPF S{x't'\pxt) for the tangent interpolation recipe. Note that while GI only 
defines a class of PSPFs in curved space, with a recipe for propagating the momentum 
needed to pick out a single PSPF, the TI procedure defines one unique path and thus 
only one PSPF. 

In flat space, it is easy to verify that the TI prescription leads to the PSPF (pTj), 
by substituting ( pO|) into the flat-space equivalent of (^), 

1 (x' - x)'^ 

S{x't'\xpxt) = p- {x - x) + -—^ -. (22) 



The first correction term leads to the form (|TO|,|TT]) with the value A = — 1. This means, 
from the results of the previous section, that it satisfies lOF, at least to the first non- 
trivial order in the distances involved. 

In fact, the TI PSPF satisfies the involution property to all orders, because the TI 
prescription obeys lOP, even over finite intervals. lOP can be verified for fiat space as 
follows (Fig. IpD): Since the sum of the principal functions is given by ([TsD , the values 
of x' and p' which extremize it are again given by (|19|). The auxiliary point x' for the 
second interval is given, analogously to (pOD, by x' = x" — p'(t' — t) = x'; since x' and 
x' coincide, there is no discontinuity in the path from x' to x", which is just a classical 
phase space path. The point x = x' — p(t' — t) = x" — p{t" — t) to which the initial jump 
at time t is made is in the same location whether the final condition is taken to be x" 
at t" or x' [given by (|19D] at t'. 



4. CONCLUSIONS AND OUTLOOK 



We have shown that the involution property described in Sec. ^ is obeyed by the phase 
space skeletonization procedure which we defined in Sec. ^ and named tangent inter- 
polation. According to prior results, summarized in Sec. |1|, using this skeletonization 
rule to construct a phase space path integral will induce a particular measure on the 
corresponding configuration space path integral. In this way, involution singles out 
a preferred realization of the skeletonized path integral. It is one which leads to a 
Schrodinger equation whose Hamiltonian differs from the simplest form, proportional 
to the covariant Laplacian, by a term proportional to the scalar curvature of the 
spatial manifold. Understanding the implications of this result rests on two as yet open 
questions. First, is there any physical explanarion for the modified Hamiltonian, and 
second, what is the significance of the involution property for the path integral prescrip- 
tion? Involution of the Hamilton principal function enables composition of propagators 
constructed via the configuration space path integral, for any choice of measure. Per- 
haps involution of the phase space principal function is related to composition of path 
integrals defining momentum space propagators. Further investigation is called for. 
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